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_^ ABSTRACT 

> 

^ Employing our previous framework to treat non-perturbative effects self- 

^ consistently, including duality violations, we update the determination of the 

cn strong coupling, a^, using a modified version of the 1998 OPAL data, updated 

^r^ to reflect current values of exclusive mode hadronic r decay branching fractions. 

O Our best n/ = 3 values from the updated OPAL data are asijnl) = 0.325 ±0.018 

and as(m^) = 0.347 ± 0.025 in fixed-order and contour-improved perturbation 
theory, respectively. 

To account for non-perturbative effects, non-linear, multi-parameter fits are 
necessary. We have, therefore, investigated the posterior probability distribution 
Cd of the model parameters underlying our fits in more detail. We find that OPAL 

data alone provide only weak constraints on some of the parameters needed 
to model duality violations, especially in the case of fits involving axial vector 
channel data, making additional prior assumptions on the expected size of these 
parameters necessary at present. We provide evidence that this situation could 
be greatly improved if hadronic spectral functions based on the high-statistics 
BaBar and Belle data were to be made available. 
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I. INTRODUCTION 

In a previous article [1], henceforth referred to as PI, we developed a new framework 
for the determination of the strong coupling, as{m'^), from non-strange vector {V) and 
axial-vector (A) hadronic r-decay data. The new framework starts from the usual finite- 
energy sum-rule (FESR) analysis, but improves this approach in two ways with regard to 
the small, but quantitatively significant non-perturbative corrections present in the theoret- 
ical representation of the FESR spectral integrals below the r mass. First, contributions 
from higher orders in the operator product expansion (OPE) are taken into account self- 
consistently. Second, in view of the fact that duality violations (DVs) are clearly present in 
the experimental spectral distributions, we use an explicit parametrization of violations of 
quark-hadron duality in our fits. As explained in detail in PI, these two improvements are 
intricately connected: estimates of the non-perturbative contribution to the sum rules with 
controlled errors cannot be obtained without taking both of these effects into account. 

Our framework was tested in VI by applying it to data from the OPAL collaboration [2]. 
We showed that fits to the data using this new framework are indeed feasible in practice. 
The resulting value for as acquires larger errors than seen in previous extractions of ag 
from hadronic r decays. The most important reason for this is that, in order to take DVs 
into account, our fits necessarily contain more parameters, while we are limited to presently 
available data. 

More recent data are in principle available. First, there are the ALEPH data [3], up- 
dated in 2005/08 [4, 5]. Presently, use of the 2005/08 ALEPH data is questionable because 
correlations due to unfolding were inadvertently omitted in the 2005/08 ALEPH update 
and hence from the publicly available covariance matrices [6]. Alternatively, more precise 
spectral functions can in principle be extracted from BaBar or Belle data. This would be 
very interesting, because one expects such spectral functions to have significantly smaller 
errors in the energy region near the r mass important for the extraction of ag- We will argue 
in this article that it should be possible to determine the non-perturbative contributions to 
the sum rules, and thus a^, with much smaller errors were such data to become available. 

It is nevertheless possible to make some progress with the OPAL data beyond the results 
presented in PI. The reason is that in Ref. [2] the normalizations of the exclusive r decay 
modes, as well as the values of a number of physical constants (such as the r mass, the 
electronic branching fraction B^, etc.) were taken from the 1998 Particle Data Group (PDG) 
tables. More precise values for these branching fractions and constants are now available 
from Refs. [7, 8], and, using these, it is thus possible to, at least partially, update the OPAL 
spectral functions. Carrying out this update, and refitting the resulting modified weighted 
spectral integrals using the methods developed and tested in VI is the primary aim of the 
present article. 

We have also investigated the probability distribution of the model parameters that un- 
derlies the various fits to the OPAL data in much more detail, using a Markov-chain Monte 
Carlo (McMC) code in order to map out the a posteriori distribution. This is useful since 
the fits we perform are non-linear in the parameters, so that not much is known a priori 
about the shape of the probability distribution. This exploration helps with understanding 
various potential instabilities in the fits (as already detected in VI), as we will discuss in 
detail below. 

As in VI, we carry out the analysis using both fixed-order perturbation theory (FOPT) 



and contour-improved perturbation theory (CIPT) [9].^ In both cases, we find that the 
central values for Ug increase compared to the values found in PI, though the two sets of 
values are consistent within errors. The errors themselves stay approximately the same, 
which is no surprise, because they are primarily determined by the errors on the OPAL 
spectral data. 

In Sec. II we briefiy review the essentials of the theory needed to understand the para- 
metrization used in our fits to the OPAL data. In Sec. Ill we explain in detail how we used 
recent results from the Heavy Flavor Averaging Group (HFAG) [7] to update the OPAL 
spectral functions. Some details are relegated to an appendix. Then, in Sec. IV we discuss 
what can be learned from the posterior probability distribution obtained with the McMC 
code. We present the results of our fits in Sec. V and summarize them in Sec. VI. In Sec. VII, 
we argue that the reduction of errors on the spectral functions expected from the BaBar or 
Belle data are likely to be of significant help in reducing the non-perturbative uncertainties. 
Section VIII contains our conclusions. 



II. THEORETICAL PARAMETRIZATION 

We start with a very brief review of the theory underlying our fits, referring the reader 
to PI for more details. Our fits are based on FESRs of the form [15, 16] 

^-'-' - [ 'i -w ^- '^' - -h L ^ '"''' "- '^* ■ <"' 

where the weight w{s) is a polynomial in s, and Ily,^ (s) with s = q'^ = —Q^ is defined by 

d'xe^'^^ (0|T {J,{x)Jm} |0) = {q,qu - q'g.u) U^'^'\s) + q'g,.U^'\s) . (2.2) 

Here J^ is one of the non-strange V or A currents U'~f^d or u'~f^'~fr,d, and the superscripts (0) 
and (1) label spin. 

The spectral functions Py/^ are taken from OPAL [2], and the integral on the left-hand 
side of Eq. (2.1) is then approximated by a sum over bins, with sq G [smin, Smax], which is 
our fitting interval. These data do not contain the pion pole, which needs to be added by 
hand. Other (pseudo-) scalar contributions are numerically negligible, being suppressed by 
two powers of the light quark masses,^ 

pf(s) = 0[{m^-m,)'], (2.3) 

pf{s) = 2fl {5{s - ml) - 5{s)) + 0[{m^ + ruaf] . 

In our fits, we will use the value /^ = 92.21 ±0.14 MeV [8]. The right-hand side of Eq. (2.1) 
provides the connection to theory, and is parametrized in terms of the strong coupling 



^ For recent investigations of these two resummation schemes, see Refs. [10-14]. 

^ The second (5-function in ^^4(5) comes from the kinematical singularity in Eq. (2.2). However, the 
combination p^^+^> is free from kinematical singularities. 



as{iTi^), the OPE condensates, and a parametrization of the DV part of Ily,^ (s). We write 
(for both V and A) 

n(^-'') is) = 4^,°) (.) + ng;,°) (.) + 4 V"°) (.) , (2.4) 

with the subscripts "pert," "OPE," and "DV" denoting the perturbative, OPE (of dimension 
larger than zero), and DV contributions to 11*^^"'"°^ (s). 

The perturbative part of the right-hand side of Eq. (2.1) can, by partial integration, be 
written in terms of the perturbative Adler function 

D'it'\s) = _,<p!l (2.5) 



ds 



n+l / \ fe-l 



n=0 k=l ^ ^ 

where ag^fi"^) = a;s(/i^)/7r. Since D{s) is independent of /i, we can choose (for instance) 
fi^ = So in Eq. (2.1), which corresponds to the FOPT scheme, or /i^ = — s, which corresponds 
to the CIPT scheme [9]. We will use values for the coefficients c„i calculated in Ref. [17] up 
to order n = 3 and in Ref. [18] up to order n = 4; for C51 we use the estimate C51 = 283 ± 283 
of Ref. [10]. The values of Cnk for k > 1 follow from the Cni using a renormalization-group 
analysis based on the fact that the Adler function is independent of /i [19]. 

The (higher-dimension) OPE contribution can be expressed in terms of the OPE coeffi- 
cients CD=2k as 

nS^'W = f:M- (2-6) 

k=l ^ ' 

In our fits we will set C2 = (it is purely perturbative and suppressed by two powers of the 
light quark masses),^ and we will treat C4, Cq and Cs as constant, neglecting logarithmic 
s dependence; we will have no need for the coefficients Cd>8- To leading order in as, and 
ignoring tiny isospin-breaking effects and perturbative light-quark mass contributions, C4 is 
the same in the V and A channels; this is not the case for Cq and Cg. For a more detailed 
discussion, including references, see VI. 

Finally, the DV contribution to the right-hand side of Eq. (2.1) can be expressed in terms 
of the DV part of the spectral function 



p^^{s) = -lmU^^^'\s), (2.7) 



as [22] 



n 



V^so) = ~<f - w{s) Ui;^'\s) = - r -«;(s)p°V(s) . (2.^ 



27ri Jis\=so So Jso ^0 

In a slight variation on Ref. [22], we parametrize PyXi ^^ 

PvyAis) = exp {-5v/A - 1v/as) sin {ay/A + (3v/as) ■ (2.9) 

^ For a alternative view of the D = 2 contribution in this context, see Refs. [20, 21]. 



This adds four new parameters per channel, in addition to Us and the OPE coefficients, 
to the fits to Eq. (2.1). The interval [smimSmax] has to be chosen such that the expres- 
sions (2.5), (2.6) and (2.8) with (2.9) provide an accurate representation of the right-hand 
side of Eq. (2.1) over the whole interval. The ansatz (2.9) was developed in Refs. [23, 24], 
based on the earlier ideas of Ref. [25].^ 

In Eq. (2.9), we have traded the parameters kv/a of VI for the parameters Sy/A] they are 
related (for both V and A) by 

K = e-^ . (2.10) 

The reason for making this change is that the fit errors on 6 are much more symmetric 
than those on k. The (strong) correlations between k, and 7 in each channel correspond to 
correlations between 6 and 7 which are much closer to linear. 
In this article, as in VI, we will employ the weights 

M^) = 1 , (2.11) 

W2{x) = I — x'^ , 

w^ix) = (1 - a;)2(l + 2a;) , 

X = s/sq ■ 

The weight W3 corresponds to the (spin-1) kinematic weight that appears in the hadronic 
branching ratio Rr- Note that 

4:t.(^o) = l27r'SEw\Vu,\'lif:l,{so) (2.12) 

is, for So = iTT-Ti equal to the (1+0) contribution to the ratio of the non-strange hadronic decay 
width and the electronic decay width of the r. In the following, we will find it convenient to 
distinguish between Jcx (so); denoting the left-hand side, and J,:^ (sq); denoting the right- 
hand side of Eq. (2.1). 

If we choose C4, C^y/A and Csy/A constant, it follows that none of these coefficients 
contribute to I^^° , only Cgy/^ contribute to /th ^ , and both C^y/A and Cgy/A contribute 
to /th • Since we will not use weights of degree larger than 3, there is no need to consider 
the OPE coefficients Cd with D > 8. 

Weights w{x) which are functions of the dimensionless variable x = s/sq are chosen in 
order to facilitate the separation of OPE contributions to I^^ (sq) having different D = 2k 
which, with this choice, scale as l/sg. While non-perturbative contributions are small at 
the scales of typical r-decay analyses, at the level of precision claimed in recent as determi- 
nations they are definitely not negligible. For example, almost the entire difference between 
the results of Refs. [5] and [27] are due to differences in the fitted non-perturbative contribu- 
tions. As discussed in detail in Ref. [27] and VI, taking advantage of the So-dependence of 
the moments lex (sq) is crucial for properly constraining such higher-D contributions. For 
further discussion of the selection of the particular set of weights chosen above we refer the 
reader to PI. 



^ The parametrization of DVs is also discussed in Ref. [26] . 

^ We have checked the influence of higher-order as corrections to the D ^ 4 contributions in the OPE. 
Numerically, the differences are tiny, and can safely be neglected. For more discussion, see VI. 



III. THE OPAL DATA UPDATE 

The 1998 OPAL inclusive pv and pa distributions were constructed as sums over exclu- 
sive mode distributions. In this process, the distributions of the three main hadronic modes 
in each channel (tt^tt'^, tt^Stt^ and 7r^7r+7r^7r° for the V channel and 7r^27r°, tt^tt+tt^ and 
7r~7r+7r~27r'^ for the A channel) were explicitly measured, while the small residual contribu- 
tions associated with other modes were typically Monte Carlo generated using TAUOLA 2.4 
[2]. The normalizations of the exclusive modes (residual or not) were, however, not mea- 
sured by OPAL, but rather fixed by the 1998 PDG values for the exclusive-mode branching 
fractions. Significant improvements to these branching fractions have been made since 1998. 

Since the distributions for the main exclusive modes noted above are publicly available, 
it is possible to update the dominant contributions to the inclusive V and A distributions by 
simply rescaling these contributions with the ratio of the new and old branching fractions for 
these exclusive modes. Unfortunately, this is not the case for the residual mode contributions 
since the individual Monte-Carlo-generated residual exclusive-mode distributions are not 
publicly available. The distribution for the sum of residual modes in each channel is, however, 
reconstructable from the publicly accessible inclusive and exclusive mode distributions. This 
distribution may then be updated in an averaged sense by computing the new and old 
versions of the sum of residual-mode branching fractions and rescaling the old combined 
residual-mode distribution by the ratio of these results. Since different exclusive modes have 
different s-dependent distributions, this average updating of the residual distributions is not 
perfect. Fortunately, however, the residual modes do not play a major role in the spectral 
functions in the kinematically accessible region (accounting, for example, for only 2.6% of the 
inclusive branching fraction in the V channel and only 1.7% in the A channel). The average 
rescaling required for the combined K-channel residual branching fraction turns out to be 
small (reducing the OPAL combined residual-mode branching fraction sum by only 1.7%). In 
contrast, the HFAG version of the A-channel combined residual branching fraction is 1.394 
times the corresponding OPAL value, making the average residual-distribution rescaling 
procedure much safer for the V channel than it is for the A channel. 

We perform the updates of both exclusive mode distributions and the combined residual 
mode distributions using branching fractions from a recent unitarity-constrained HFAG fit.^ 
The particular fit we employ is that incorporating Standard Model expectations based on 
7r^2 and K^2 data for B[t ^ nUr] and B[t ^ Kvr] in addition to the results for these 
branching fractions measured directly in r decays [7].^ 

It is important to note that the conventions for quoting the various exclusive branch- 
ing fractions are not identical for OPAL and HFAG. HFAG quotes cjvr", wk^-k^ and rj-K^-K^ 
branching fractions corresponding to all u and rj decay modes, and excludes uj and rj substate 
contributions in quoting branching fractions for all other modes. In contrast, for OPAL, (i) 
the quoted 7r~7r+7r~, 7r~7r"'"7r~7r° and 7r~7r+7r~27r° branching fractions include, respectively, 
WTT", WTT" and W7r~7r'^, and wtt^tt^ and rjii^Ti^ components, and (ii) the U7r~ and wtt^tt^ 
branching fractions are quoted excluding cu — )> Svr contributions. With these conventions, 
the tabulated OPAL exclusive branching fractions and distributions include small '"wrong- 
current contaminations" associated with isospin-breaking u — )■ tt^tt" and r] — )■ 7r"'"7r~7r° 



^ The updated OPAL data are available on request. 

^ We refer to littp://www. slac.stanford.edu/xorg/hfag/tau/lifag-data/tau/2009/TauFit_Mar2011 
/BB_PiKUniv/ConstraincdFit.pdf for details. 



decays. Explicitly, oo — )> tt+tt" decays cause the K-current-induced ujtx~ mode to populate 
the nominally A-current tt^tt+tt'^ distribution and the A-current-induced wk^ti^ mode to 
populate the nominally V-current Tr^Tr+Tr^Tr" distribution, while rj — )■ tt+tt^tt'^ decays cause 
the l^-current-induced rjii^Ti^ mode to populate the nominally A-current 7r^7r+7r^27r° dis- 
tribution. In forming the inclusive V and A spectra, OPAL corrects for this contamination 
by including an appropriate negatively weighted version of the relevant non-o; — )■ Stt and 
non-r] — )■ Svr utt~ , U7r~7i^ and rjir'n^ distributions in the wrong-current inclusive distribution 
sum. The negative weights employed by OPAL were determined using the 1998 PDG values 
for the branching fractions of the relevant rj and u decay modes. 

In order to perform the rescaling of the OPAL exclusive-mode distributions, the relevant 
updated wrong-current contaminations must be added to the HFAG exclusive branching 
fractions. The HFAG-updated un' (excluding u — )■ Stt), un~7i^ (excluding u — )■ Stt), and 
ri7T~TT^ (excluding rj — )■ Sir) branching fractions, and updated negative- weight, wrong-current 
contamination corrections must, analogously, be incorporated in the updated version of the 
combined residual mode branching fractions in both channels. These updates are performed 
using the HFAG exclusive branching fractions, together with 2010 PDG results for the 
relevant 77 and u branching fractions. Numerical details may be found in the Appendix. 

OPAL has also tabulated the covariance matrices for the three main exclusive modes in 
each channel, as well as the VV, VA and AA covariances for the inclusive V and A distri- 
butions. The absence of information on the covariances among the different exclusive mode 
distributions limits our ability to update the inclusive VV, VA and AA covariances. Updates 
for improvements in factors such as B^. and V^d which enter when converting the differential 
branching fraction distributions, dBy / a{s) / ds , to the corresponding spectral functions, can, 
however, be performed. Details on carrying out this procedure may also be found in the 
Appendix. 



IV. THE POSTERIOR PROBABILITY DISTRIBUTION 

The fit functions used in the sum rules (2.1) are non-linear in Q;s(m^) and the DV param- 
eters. It is therefore not obvious what the posterior probability distribution of the model 
parameters looks like, even if we assume the data errors to follow a (multivariate) gaussian 
distribution. 

In order to study this distribution, we have used an McMC code, Hrothgar [28], in order 
to generate the conditional probability distribution, which we take to be proportional to 
exp[— x^(p)/2], given the data, where p represents the array of fit parameters. With the 
data fixed, these parameters are varied stochastically, and a Metropolis-Hastings accept- 
reject step is used to generate a statistical picture of the probability distribution. In this 
section, we will describe our findings in more detail for the case of a fit to the FESR with 
weight Wq, first in the V channel. 

The McMC code generates points in the 5-dimensional space spanned by the five pa- 
rameters asiml), (5y, 7y, ay and /Sy, and also computes the value of x^ at each of the 
generated points. These points are distributed following exp[— x^(p)/2], with y^ip) evalu- 
ated on the (updated) OPAL data (including the full covariance matrix) and the values of 
the parameters p at these points. 

The probability distribution thus obtained can be projected onto two-dimensional planes. 
In Fig. 1, we show x^ as a function of asijnl), choosing Smin = 1-5 GeV^, using FOPT for 
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FIG. 1: as(m^) versus x^ J V channel, with Smin = 1-5 GeV^ (200,000 points). 




FIG. 2: Two-dimensional contour plots showing as{m^) versus Sy and 7y versus 5v ■ Left panel: 
projection onto as{m'^) — 6v plane. Right panel: projection onto jy ~ ^v plane. V channel, 
Smin = 1-5 GeV . Blue (darker) areas and green (lighter) areas contain 68%, respectively, 95% of 
the distribution. 



the perturbative part.^ Since for each as{ml) points with many different values for the other 
four parameters are generated stochastically, the distribution appears as the cloud shown in 
the figure.^ 



^ The distribution for CIPT looks essentially the same, except that the projections shown in Figs. 1 and 2, 

left panel, are shifted to the right by an amount ~ 0.02. 
^ If a new point is rejected by the accept-reject step, the old point is retained. Therefore each point in the 
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FIG. 3: x^ versus as{'m'^), thinning out the integrated data by factors 1 (upper left), 2 (upper 



right), 3 (lower left) and 4 (lower right); 
(200,000 points). 



V and A channel combined, with Sr, 



1.5 GeV2 



Figure 1 shows a bi-modal distribution, with one local minimum near as{m^) = 0.28, 
and a global minimum near ^^(m^) = 0.32, with a difference in the locally minimal values 
of x^ equal to about 1.6. As a consequence, a standard x^ minimization, which estimates 
the parameter covariance matrix from the hessian at the (global) minimum, ^'^ will miss the 
other local minimum entirely. 

The origin of the problem appears to be the fact that Sy is not well constrained by the 
data. This can be seen in Fig. 2, which shows the projections onto the as{'m'^) — 6 and 
7 — (5 planes, in this case as contour plots showing the regions containing 68% (blue) and 
95% (green) of the distribution. The right panel shows a very strong correlation between 
the two parameters, 6v and 7v, which together control the "strength" of the DV part of the 
spectral functions in the low-s part of our fitting windows, cf. Eq. (2.9). Clearly, external 
input is required to narrow down which part of the distribution is most likely to correspond 
to physics. This will be discussed in Sec. V. 

Figure 3 shows analogous results for a V and A channel combined (V^A) fit, again for 
Smin = 1-5 GeV^, again using FOPT. In this case, there are 9 fit parameters, and the figure 
shows a projection of the 10-dimensional space spanned by the 9 parameters and x^- The 
four panels show a fit of the V^A FESR with moment wo, similar to the V channel fit shown 



plot may represent multiple points in the ensemble. 
10 Qj. fj-oni the minimum value of x^ plus one. 



in Fig. 1. In the upper left panel we used the full set of sq values corresponding to the right 
end-points of all bins starting from Smin = 1-504 GeV^, whereas the other three panels show 
fits in which the Sq values employed have been thinned out by a "thinning factor" n, chosen 
equal to 2, 3 and 4, respectively, in the upper right, lower left and lower right panels. ^^ 
Contour plots for the combinations ^^(m^) — 6y and 7y — Sy look very similar to those 
shown for the V case in Fig. 2. 

Again, as in Fig. 1, there appear to be two local minima, one centered around ^^(m^) = 
0.315, and one centered around ^^(m^) = 0.28. However, in this case the two minima are 
much closer to being degenerate than in the V-channel fit. For n = 4 it is difficult to discern 
two separate minima; the two minima appear to merge. 

Similar behavior as a function of the thinning factor is observed in fits to only the V 
channel as well, but the two minima are always more clearly separated, as in Fig. 1. This 
may explain why the VSzA combined fits in VI were found to be less stable than V channel 
fits. Figure 3 may also explain why fits with n = 3 led to more stable results in VI, since for 
n = 3 the two minima appear to be somewhat more clearly separated than for other values 
of n}'^ 

One should also bear in mind, in assessing the relative reliability of the results of the V 
channel and combined V^A channel fits, that the much larger average rescaling of the sum- 
of-residual-modes branching fraction in the A channel makes the updating of the spectral 
function much less reliable in the A channel than it is in the V channel. 

For both the V-onlj and the V^A fits, we have also studied the behavior of the x^ 
distribution as a function of Smin- We find that by lowering Smm, the minimum at the lower 
value of asi'm'^) "moves up," i.e., the value of x^ at that local minimum increases relative 
to that at the other minimum. 



V. FITS 

In this section, we present the results of fits to a range of different moments, obtained 
by minimizing various different "fit qualities," (positive-definite quadratic forms in the dif- 
ferences between theory and data). In Sec. VA we discuss fits to the FESR with moment 
jiwo)^ whereas in Sec. VB we will consider simultaneous fits to FESRs with the moments 
j(^o,2,3)^ using the weights of Eq. (2.11). 

In the first case, we choose the fit quality to be the standard x^, already discussed in the 
previous section, constructed with the complete (updated) covariance matrix. In the second 
case, in which we combine more than one moment, it turns out, as discussed in VI, that 
the correlations are too strong to allow for a fit based on the standard x^ function. ^^ We 
therefore employ a somewhat simpler fit quality Q^. Working with a set of values of sq, {sq} 



^^ If the thinning factor is equal to n, we use every n-th value of lex" (sq) in the fit, see also VI. We 

emphasize that all data are used: only integrated data are thinned out. 
^^ We have checked that the behavior of the posterior probability distributions with the non-updated data 

we used in 7^1 is very similar to what we find with the updated data. 
^^ When more than one weight is employed, the correlation matrix for the full set of weighted spectral 

integrals, labelled by the weights and sq values employed, acquires a number of machine-precision zero 

eigenvalues. 



10 



in some fitting window, we define 



with Cw the covariance matrix for moments with fixed weight w and sq running over the 
chosen fit window. ^^ The fit quahty Q^ is thus similar to a standard x^^ but with cross- 
correlations between different moments omitted. Treating Q^ as if it were the standard x^ 
would thus lead to incorrect errors and covariances for the fit parameters. To take the cross- 
correlations properly into account, errors and covariances for Q^-based fits are obtained 
using the linear fluctuation analysis described in the Appendix of PI. 

In view of (i) the results of Sec. IV and (ii) the fact that the updating scheme is much 
more reliable for the V channel than for the A channel OPAL data, we will use the l^-channel 
fits of Sec. VA to determine our central value for asirri^). The remaining fits are used only 
to investigate whether our fit function, which parameterizes DVs using Eq. (2.9), provides 
a good description of the data for the moments W2 and W3 as well. The issue of the choice 
of weights is discussed in more detail in VI. 

As we have seen in Sec. IV, the posterior probability distribution generally has a rather 
complicated structure, showing almost always two fairly close but different minima. We 
thus need to address the question which minimum is more likely to correspond to a physical 
solution. The situation is more complicated in the case of V&:A fits, for which the two 
minima are essentially degenerate. We will argue that the minimum corresponding to the 
larger value of Q;s(m^) is more likely to correspond to the correct physics. 

First, there is evidence for this choice from the fits themselves. We note that, in the 
V-onlj case, the minimum of x^ ior the larger value of as{'m'^) is always the lower one by a 
significant amount, cf. Fig. 1. This is confirmed by fits with lower values of Smin, for which 
this separation becomes more pronounced. 

We may also refer to the model study of Ref. [24] , which led to the form of the ansatz (2.9) 
used to parameterize DVs. It was shown there that the model underlying this ansatz leads 
naturally to the following values for the parameters: 

(F'^\ 11 

— ~4 and 7 ~ — — ~ 0.3 GeV"^ , (5.2) 

A2y ' iVcA2 ^ ' 

where F ~ 0.1 GeV is a typical value for a resonance decay constant, and A ~ 1 GeV is a 
typical QCD scale. Figure 2 shows that for such values the global x^ minimum, which occurs 
at the larger value of 5y, and thus at the larger value of q;s(?72^), is preferred. (We will see 
below that the estimates of Eq. (5.2) are less well satisfied for the A channel.) Henceforth, 
we will refer to the minima at larger values of 5v as "physical" minima, and to those at 
smaller values of 5v as "unphysical." 



A. Fits with wq 

In Table 1 we show V channel fits of l[t'\sQ) to /ir^(so) (c/ Eq. (2.1)), for sq G 



The fit quality Q^ corresponds to Qli^^i, defined in T'l. 
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Smin 


dof 


xVdof 


as 


5v 


IV 


av 


Pv 


1.3 
1.4 
1.5 
1.6 
1.7 


53 
50 
47 
44 
41 


0.41 
0.33 
0.34 
0.35 
0.35 


0.338(18) 
0.326(16) 
0.323(16) 
0.325(18) 
0.323(19) 


3.91(62) 
4.11(63) 
4.21(62) 
4.04(86) 
4.37(99) 


0.27(43) 
0.16(41) 
0.12(40) 
0.20(53) 
0.05(53) 


0.53(54) 
-0.29(68) 
-0.48(79) 
-0.37(87) 
-0.48(91) 


2.89(29) 
3.29(35) 
3.38(40) 
3.33(44) 
3.37(44) 


1.3 
1.4 
1.5 
1.6 
1.7 


53 
50 
47 
44 
41 


0.43 
0.34 
0.34 
0.36 
0.37 


0.360(32) 
0.349(25) 
0.345(24) 
0.347(26) 
0.344(25) 


3.47(64) 
3.84(65) 
3.99(64) 
3.82(90) 
4.2(1.1) 


0.53(47) 
0.30(44) 
0.22(42) 
0.32(56) 
0.12(57) 


0.57(58) 
-0.31(67) 
-0.54(77) 
-0.45(86) 
-0.57(90) 


2.83(32) 
3.28(35) 
3.39(40) 
3.35(44) 
3.40(44) 



TABLE 1: Standard x^ fits to Eq. (2.1) with w{s) = 1, V channel. FOPT results are shown 
above the double horizontal line, CIPT results below. Errors are standard x^ errors; ^y and fiy in 
GeV-2. 



[smini Smax\) with Smax = 3.136 GeV^ and varying Smin-^^ In all the fits contained in this 
table, we have used initial parameter estimates which roughly correspond to the physical 
minima, i.e., the minima corresponding to larger values of 6v found with the McMC code. 
There is excellent stability for Smin ranging from 1.4 to 1.7 GeV^. In a slight deviation 
from VI, we will use the average of the fits with Smin = 1-4, 1.5 and 1.6 GeV^ to determine 
as{fn'^). Since the fit is non-linear, one expects the fit errors to be asymmetric. For instance. 



for the FOPT fit with s. 



1.5 GeV^ we find 



as{ml) 


n QOQ +0.016 


Sv 


= 4.2ii°:i. 


Iv 


= 0.12lO:iGeV 


ay 


= -o.48;°:sl , 


Pv 


= 3.38+°:^^ GeV 



(5.3) 



We note that the error on as{m\) is nearly symmetric, and that the error on 5y is much 
closer to symmetric than the error on ny = exp{—Sy) in VI. A typical parameter correlation 
matrix, that for the FOPT fit with Smin = 1.5 GeV^, is shown in Table 2. Results for other 
values of Smin, or for CIPT fits, show the same pattern. 



Recalling our choice to obtain a central value by averaging results for s^ 
1.6 GeV^, we obtain from these fits for a, at the r mass the results 



1.4, 1.5 and 



asiml) 



0.325 ± 0.016 ± 0.002 ± 0.007 (FOPT) , 
0.347 ± 0.024 ± 0.002 ± 0.005 (CIPT) . 



(5.4) 



^^ This value of Smax corresponds to the highest bin available in the OPAL data; the bin width is 0.032 GeV^. 
In the axial channel the highest bin available corresponds to s^ax — 3.104 GeV^. 
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Ois ^V IV Civ Pv 


as 


1 0.68 -0.67 0.74 -0.68 


Sv 


0.68 1 -0.99 0.47 -0.44 


IV 


-0.67 -0.99 1 -0.49 0.45 


ay 


0.74 0.47 -0.49 1 -0.98 


Pv 


-0.68 -0.44 0.45 -0.98 1 



TABLE 2: Parameter correlation matrix for the FOPT fit with Sr 



1.5 GeV^ shown in Table 1. 




1.5 2.0 2.5 

^0 (GeV2) 



0.0 0.5 1.0 1.5 2.0 2.5 3.0 

^0 (GeV2) 



FIG. 4: Left panel: comparison of lex'^ (sq) and I^i^'^ (sq) for the Smin = 1-5 GeV V- channel fits 
of Table 1. Right panel: comparison of the theoretical spectral function resulting from this fit with 
the experimental results. CIPT fits are shown in red (dashed) and FOPT in blue (solid). The 
(much flatter) black curves represent the OPE parts of the fits. The vertical dashed line indicates 
the location of Smin- 



The first error is the Sr 



1.5 GeV^ fit error shown in Table 1, the second the variation 



of the central values over the Smin = 1-4 — ?■ 1.6 GeV^ averaging window, and the third the 
result of the ±283 variation of C51 about its central value C51 = 283. 

Figure 4 shows the i&o-FESR fit for Smin = 1-5 GeV^ (left panel), and the corresponding 
theoretical curves for the spectral function in comparison with the (updated) experimental 
data (right panel). Agreement with data is good in the full fit window sq > Smin = 1-5 GeV^. 
The black curves show the OPE parts of the theoretical curves, i.e., the curves obtained by 
removing the DV contributions from the blue and red curves. Clearly, DVs are needed to 
give a good description of the data for Iex° and the spectral function itself. We emphasize 
that the right panel of Fig. 4 is not a fit; only the moments I^^^^Sq) were used in the fits 
reported in Table 1. 

Fits with weight wq to the combined V and A channels are tabulated in Table 3, where 
again initial parameter estimates were chosen close to the physical minima. In view of our 
findings of Sec. IV for this case, we chose the thinning factor n equal to 2. We do not 
show plots of these fits, or of the corresponding spectral functions, because they look very 
similar to those shown in Fig. 4 and the corresponding figures in VI. The CIPT fit with 
Smin = 1.7 GeV^ appears to correspond to an unphysical solution of the type discussed in 
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Sec. IV; we did not find a physical solution in this case. 

Following the same prescription as for Eq. (5.4), we obtain for Ug the values 



as 



m. 



m^ 



0.319 ± 0.015 ± 0.007 ± 0.005 (FOPT) , 
0.338 ±0.021 ±0.010 ±0.004 (CIPT) 



(5.5) 



from the VS^A fits. The errors have the same meaning as in Eq. (5.4). The values in Eq. (5.5) 
are in good agreement with those of Eq. (5.4). It should, however, be kept in mind that (i) 
the physical and unphysical minima of the x^ function are close to degenerate for these fits 
(c/. Sec. IV) and (ii) the averaged rescaling of the sum-of-residual-modes part of the OPAL 
spectral functions is considerably less reliable for the A channel than for the V channel. 



Smin 


dof 


xVdof 


as 


W/A 


IV/A 


Oiy/A 


f3v/A 


1.3 


49 


0.58 


0.327(12) 


3.71(55) 
1.62(86) 


0.38(39) 
1.66(55) 


0.24(51) 

2.48(77) 


3.00(30) 
3.60(46) 


1.4 


46 


0.47 


0.325(11) 


4.28(44) 
1.6(1.0) 


0.02(30) 
1.68(64) 


-0.54(58) 
1.8(1.2) 


3.43(31) 
4.00(68) 


1.5 


43 


0.53 


0.312(15) 


3.90(71) 

1.82(72) 


0.29(46) 
1.46(44) 


-1.02(85) 
-2.5(1.3) 


3.64(45) 

2.91(74) 


1.6 


40 


0.40 


0.320(13) 


4.23(57) 
1.56(94) 


0.05(35) 
1.64(53) 


-0.79(70) 

2.8(1.7) 


3.55(37) 
3.47(91) 


1.7 


37 


0.54 


0.312(17) 


3.7(1.4) 
0.3(1.8) 


0.42(78) 
2.15(86) 


-0.9(1.0) 

-1.7(2.1) 


3.60(51) 
2.5(1.1) 


1.3 


49 


0.61 


0.348(18) 


3.38(51) 
1.95(78) 


0.58(38) 
1.48(50) 


0.30(54) 
2.51(83) 


2.93(32) 
3.61(49) 


1.4 


46 


0.49 


0.347(15) 


4.03(46) 

1.97(82) 


0.14(32) 
1.48(53) 


-0.60(54) 
2.0(1.1) 


3.44(30) 
3.93(64) 


1.5 


43 


0.54 


0.328(21) 


3.69(77) 
1.94(71) 


0.39(51) 
1.39(42) 


-1.08(83) 
-2.4(1.3) 


3.66(45) 
2.90(72) 


1.6 


42 


0.40 


0.339(17) 


4.09(61) 
1.73(94) 


0.12(38) 
1.54(52) 


-0.90(67) 
2.9(1.5) 


3.59(35) 
3.42(83) 


1.7 


37 


0.47 


0.294(17) 


-0.9(2.8) 
0.8(1.7) 


3.3(1.6) 
1.76(78) 


3.1(2.3) 
-0.9(1.7) 


8.2(1.2) 
2.12(87) 



TABLE 3: Standard x^ fits to Eq. (2.1) for w{s) = 1, combined V and A channels. FOPT results 
are shown above the double horizontal line, CIPT results below. The first line for each Smin gives 
the V DV parameters; the second line the A ones. Every second value of sq in the range above and 
starting at Smin is included in the fits. Errors are standard x^ errors; jy/A o-nd l^v/A ^^ GeV~^ 



B. Fits with Wo, 2, 3 

In this section, we report on simultaneous fits to moments with weights Wq, W2 and Wa, 
using the fit quality Q^ of Eq. (5.1). FOPT and CIPT results are shown for the V channel 
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Smin 


dof 


QVdof 


as 


Sv 


7v 


av 


f^v 


W^Cey 


lO^Csy 


1.3 
1.4 
1.5 
1.6 
1.7 


167 
158 
149 
140 
131 


0.42 
0.33 
0.33 
0.33 
0.33 


0.307(13) 
0.313(13) 
0.315(14) 
0.317(17) 
0.318(19) 


2.68(75) 
3.37(69) 
3.74(60) 
3.42(76) 
4.26(73) 


1.10(51) 
0.63(46) 
0.40(39) 
0.59(48) 
0.14(39) 


0.29(80) 
-0.69(87) 
-0.9(1.0) 
-0.7(1.4) 
-0.8(1.3) 


2.89(47) 
3.45(49) 
3.55(57) 
3.48(75) 
3.53(68) 


-0.51(35) 

-0.47(28) 
-0.45(28) 
-0.42(39) 
-0.46(38) 


0.68(66) 
0.79(44) 
0.80(45) 
0.72(68) 
0.86(61) 


1.3 
1.4 
1.5 
1.6 
1.7 


167 
158 
149 
140 
131 


0.38 
0.30 
0.30 
0.31 
0.31 


0.362(45) 
0.349(30) 
0.345(30) 
0.347(42) 
0.344(40) 


3.55(80) 
3.85(66) 
3.97(61) 
3.71(71) 
4.21(74) 


0.47(57) 
0.30(44) 
0.24(40) 
0.38(47) 
0.13(42) 


0.53(98) 
-0.3(1.0) 
-0.5(1.3) 
-0.4(1.8) 
-0.6(1.8) 


2.85(52) 
3.28(55) 
3.39(66) 
3.33(92) 
3.40(88) 


-0.18(51) 
-0.40(33) 
-0.46(35) 
-0.42(55) 
-0.50(50) 


0.06(82) 
0.53(55) 
0.66(61) 
0.6(1.0) 
0.77(92) 



TABLE 4: Fits to Eq. (2.1) with weights ti)o,2,3) V channel, using fit quality (5.1). FOPT results 
are shown above the double horizontal line, CIPT fits below. 71/ and /3y in GeV~^, C^y in GeV^ 
and C^y in GeV^. 

in Table 4, for the same set of Smin values as before. To properly account for the cross- 
correlations between moments with different weights, errors and covariances are computed 
through the linear-fluctuation analysis of VI. We do not show any plots based on these fits, 
as they look very similar to those in PI. 

In this case, we have not carried out an investigation along the lines of Sec. IV. The 
reason is that we cannot compute a fully correlated posterior probability distribution, and 
the interpretation of the probability distribution associated with Q^ would be less clear. 
Our only reason for considering these multiple-moment fits is to check that DVs in higher 
moments, and in particular the moment with the kinematic weight, can be described by our 
ansatz, Eq. (2.9). We find that this is indeed the case. 

The fit results, reported in Table 4, are in good agreement with those of Sec. VA. Fol- 
lowing the same method as before, we obtain for ag the values 






0.315 ± 0.014 ± 0.002 ± 0.007 (FOPT) , 
0.347 ± 0.030 ± 0.002 ± 0.005 (CIPT) , 



(5.6) 



where again the errors have the same meaning as in Eq. (5.4). For the CIPT case we see 
that adding more moments has not improved the determination and, in fact, has somewhat 
increased the total error. In view of this observation, and the fact that the errors in Eqs. (5.6) 
and (5.4) were, in any case, produced using different minimizing functions, we stick with 
the standard x^ V^-channel fit results of Eq. (5.4) as our central ones. 

In Table 5 we show similar results for the VSzA analysis. In this case, we found the most 
stable results with the thinning factor n = 3, i.e., thinning out the moments /ex '^'^ (sq) 
by a factor three. Even so, we did not find a physical minimum for the FOPT case with 
Smin = 1-6 GeV^, as can be seen from the table. The distinction between the physical 
minima we found at lower values of Smin and the unphysical minima at higher values of Smin 
is very clear from the values of the DV parameters. ^^ In particular 6v and 7^ both differ 
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The value of x^ is always smaller at the unphysical mininium in these particular fits. 
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Smin 


dof 


QVdof 


as 


^V/A 


IV/A 


ay/A 


I^V/A 


lO^Qy/A 


lO^Q^M 


1.3 


104 


0.64 


0.309(9) 


2.80(74) 
2.51(43) 


1.00(50) 
1.11(28) 


-0.25(72) 
3.03(71) 


3.18(43) 
3.37(41) 


-0.54(26) 
0.56(20) 


0.82(49) 

-0.58(37) 


1.4 


98 


0.49 


0.310(11) 


3.33(72) 
2.26(50) 


0.64(48) 
1.22(32) 


-1.07(82) 
-2.89(95) 


3.66(47) 
3.17(54) 


-0.55(23) 
0.50(29) 


0.93(38) 
-0.40(59) 


1.5 


92 


0.48 


0.312(11) 


3.70(62) 
2.16(78) 


0.40(41) 
1.28(43) 


-1.23(95) 
-2.9(1.2) 


3.75(53) 
3.17(68) 


-0.52(22) 
0.53(34) 


0.92(36) 
-0.45(76) 


1.6 


86 


0.38 


0.292(14) 


-0.8(3.0) 
1.8(1.1) 


3.2(1.8) 
1.38(56) 


-0.9(1.8) 
-1.7(1.4) 


7.0(1.0) 
2.55(73) 


-1.14(18) 
-0.09(62) 


2.10(32) 
0.96(1.7) 


1.7 


80 


0.44 


0.312(17) 


3.78(89) 
-0.4(2.3) 


0.36(50) 
2.5(1.0) 


-1.1(1.4) 
-1.6(3.2) 


3.69(74) 
2.5(1.7) 


-0.52(37) 

0.27(78) 


0.91(66) 
0.5(2.4) 


1.3 


104 


0.53 


0.346(18) 


3.45(58) 
1.97(69) 


0.54(41) 
1.46(44) 


0.02(70) 
2.33(76) 


3.09(40) 
3.73(44) 


-0.38(27) 
0.71(22) 


0.43(46) 
-1.02(43) 


1.4 


98 


0.43 


0.339(17) 


3.75(57) 
2.03(60) 


0.34(39) 
1.39(39) 


-0.72(83) 
2.8(1.0) 


3.49(46) 
3.48(59) 


-0.51(22) 
0.59(26) 


0.73(38) 
-0.76(55) 


1.5 


92 


0.43 


0.337(17) 


3.89(53) 

2.13(78) 


0.26(36) 
1.34(45) 


-0.95(98) 
2.8(1.4) 


3.60(53) 
3.47(74) 


-0.55(23) 
0.58(31) 


0.82(41) 
-0.74(68) 


1.6 


86 


0.44 


0.335(23) 


3.56(77) 
1.7(1.1) 


0.45(48) 
1.51(58) 


-0.9(1.4) 
3.1(1.8) 


3.57(73) 
3.31(99) 


-0.55(34) 
0.50(46) 


0.79(65) 
-0.5(1.1) 


1.7 


80 


0.42 


0.332(30) 


3.79(84) 
-0.3(2.4) 


0.33(47) 
2.14(1.0) 


-1.0(1.7) 
-2.1(3.5) 


3.61(85) 
2.7(1.9) 


-0.58(43) 
0.30(85) 


0.88(84) 
0.2(2.5) 



TABLE 5: Fits to Eq. (2.1) with weights u'o,2,3; combined V and A channels, using fit quality (5.1). 
FOPT results are shown above the double horizontal line, CIPT fits below. 7\//a c-i^d /Sy/j^ in 
GeV~^, Cey and Cq^a in GeV^ and Csy and Cs,a in GeV^. The first line for each Smin gives the 
V channel DV and OPE parameters; the second line the A channel ones. Every third value of sq 
in the range above and starting at Smin is included in the fits. 



by a large amount between physical and unphysical solutions, much as shown in the simpler 
case displayed in Fig. 2. Averaging only the FOPT fits at Smin = 1-4 and 1.5 GeV^, and 



averaging as before the CIPT fits at s^ 



1.4, 1.5 and 1.6 GeV , we find for a. 



,m: 



as 



m^ 



m„ 



0.311 ±0.011 ±0.002 ±0.007 
0.337 ±0.017 ±0.002 ±0.005 



(FOPT) , 
(CIPT) , 



(5.7) 



with errors again as in Eq. (5.4). 



VI. SUMMARY OF RESULTS 



Through a more detailed statistical study of the data than we carried out in VI, we 
showed that our fits of the OPAL data sometimes allow for different local mimina of the x^ 
function, cf. Sec. IV. These solutions are most clearly distinguished by the values of the DV 
parameters 6v and 7y, and in the introduction to Sec. V we argued that the solutions with 
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large values of Sy (of order 4) and small values of 7y (of order 0.3) should be considered as 
physical, while the other minima, which always have small values of Sy (negative, in fact), 
and large values of 7y (typically of order 3) should not be considered physical. 

As in PI, we have used both Weinberg sum rules [31] and the sum rule for the electro- 
magnetic pion mass difference [32] to test our fit results. All three sum rules are well satisfied, 
at a level of precision similar to that found in VI. We have also confirmed, again as in VI, 
that our theoretical description of Ryj^^udi^^) (-^^- (2-12)) agrees, within errors, with data 
for So down to below 1.5 GeV^. 



A. The value of OLgirn^) 

We will choose the values of the strong coupling at the r mass obtained from the V- 
channel fit of /ex' as our best values. Our reasoning for doing so is twofold. First, while 
the simultaneous fits to multiple moments are in good agreement with this simple fit, no 
standard x^ fit is possible in this case. While we believe that the error estimates based 
on linear fiuctuation analysis are reasonable, it is less clear how they should be interpreted 
than those obtained from a standard x^ analysis. ^^ Second, including also the A channel 
makes the fits more complicated, because of the larger number of parameters. This effect 
is compounded by the much larger, and hence less certainly reliable, approximate rescaling 
that must be applied to the residual distribution in the A channel. In addition, we note that 
the only feature visible in the A channel is the ai resonance, and it is not clear whether the 
ansatzwe use to parameterize the DV part of the spectral function can be expected to apply 
to this resonance, even if we assume that the ansatz works well for higher resonances in each 
channel. Finally, related to this, we note that the typical values of the DV parameters we 
find in the axial channel satisfy the expectation of Eq. (5.2) less well. 

We thus find our best values for the strong coupling from Eq. (5.4): 

asiml) = 0.325 ±0.018 (MS, nf = 3, FOPT) , (6.1) 

as{ml) = 0.347 ±0.025 (MS, nj = 3, GIFT) , 

where we added the errors in Eq. (5.4) in quadrature. 
Running these values up to the Z mass Mz [33] yields^^ 



«.(MI) = 


= 0.1191 ±0.0022 


{MS, Hf = 5, FOPT) , 


«.(MI) = 


= 0.1216 ±0.0027 


{MS, Hf = 5, GIFT) , 



(6.2) 



where we symmetrized the resulting slightly asymmetric errors. 



^^ Of course, we have found that the posterior probability distribution has a comphcated behavior, so that 
physical input is required to decide which local minimum is physical, as discussed in detail in Sees. IV 
and V. 

^^ We evolved as to the Z mass in the same way as was done in Ref. [10]. Uncertainties in the running, 
associated with the use of 4-loop truncated /3 functions, uncertainties in the charm and bottom masses, 
and the choice of the nf — 3 ^ nf — 4 and n/ = 4 — > nf — 5 matching thresholds are negligible on the 
scale of the quoted errors. 
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B. Non-perturbative results 

In PI we estimated the relative deviation of the values found for the diniension-6 conden- 
sates from those given by vacuum-saturation approximation. To this end, these condensates, 
parametrized by Cgy/A, are expressed in terms of the quantities pi and ps by 

^6,vM = ^-V(gg)^('^;;"^'^^). (6.3) 

Vacuum saturation values for C^yjA then correspond to pi = P5 = 1. Performing the 
analogous analysis for the updated OPAL spectral functions, we find (employing {qq){'m'^) = 
-(272MeV)3 [34]) 

pi=3.1±2.0, p5 = 4.4±1.4 (FOPT) , (6.4) 

pi=3.1±1.6, p5 = 4.3±1.3 (CIPT) , 

using as representative values Cey/A and as{ml) of Table 5 for Smin = 1-5 GeV^.^^ As in VI, 
the values of pi and ps are insensitive to the perturbative resummation scheme. We note 
that pi changes sign relative to the central value found in VI, but also that, given the large 
uncertainties, there is no inconsistency between our earlier fits and those presented here. 

Analyses of the strong coupling from r decays are sometimes based on the ratio, Ry+A^ 
of the total inclusive non-strange branching fraction to the electron branching fraction B^ 
[16], 

Rv+A = NcSEw\Vud\Hl + 5p + 5Np) , (6.5) 

where Sp stands for the perturbative, and 6np stands for the non-perturbative contributions 
beyond the parton model. Determining 6p, and hence as, from Ry^A of course requires 
input for 5np- In the past, shortcomings in the methods used to obtain this input have 
led to a significant underestimate of the corresponding uncertainties. Analyses (such as 
those of Refs. [2, 4]) including additional higher-degree-weight FESRs, for example, were 
forced to assume that D > S contributions could be neglected for all additional FESRs. ^° 
Reference [27] avoided this problem, but, being unable to fit all required Z^ < 8 OPE 
parameters using an sq window within which neglect of integrated DVs was self-consistent, 
was forced to rely on external input for the gluon condensate, the renormalon ambiguity of 
which makes this external input potentially problematic. As shown in VI and Ref. [29], it 
is not possible to avoid these problems without considering lower sq and FESRs for which 
integrated DVs are not negligible in the full so-fitting window employed. The framework 
presented in VI and in the present article is the first to allow for a reliable estimate of 5np 
from such an analysis, and hence to bring these systematic issues on the theory side under 
control. Expressing the D = 6, 8 OPE terms, as well as the DV contributions to 5np, in 



^^ We neglect the errors on a., and {qq). 

"^^ This assumption has since been tested (and found to be poorly satisfied) by comparing the SQ-depcndence 
of the fitted theory side to that of the experimental data side of the various additional FESRs [27]. 



18 



terms of 5^^\ 5^^^ and 5^^ , respectively, we obtain from our fits with Smin = 1-5 GeV^ 






0.0 ± 1.9) ■ 10^2 , 5^^) = (- 3.7 ± 7.6) ■ 10"^ 

0.1 ±1.0) -10"^ (FOPT) , 

0.1 ± 1.8) • 10"^ , 5^^) = (- 0.6 ± 7.6) ■ 10~^ 

0.6 ± 1.4) • 10^3 (CIPT) . 



(6.6) 



In the case of FOPT, the corresponding correlation matrix is found to be: 





5(6) 5(8) ^DV 


5(6) 


1 -0.98 0.59 


5(8) 


-0.98 1 -0.54 


jDV 


0.59 -0.54 1 



TABLE 6: Correlation matrix for the quantities S^^\ S'-^^ and (5°^ of Eq. (6.6) (FOPT). 



While individually for the V and A channels the hierarchy of the non-perturbative terms 
is such that dimension-6 is the largest, dimension-8 smaller, and the DV contribution the 
smallest, due to strong cancellations in the D = 6 and DV contributions for the sum \^ -|- A, 
the D = 8 contribution turns out to be dominant. However, in view of the large uncertainties, 
it is impossible to conclude that these cancellations will also persist once more precise data 
are available. 

Combining the OPE contributions as well as the DV term of Eq. (6.6) including correla- 
tions, the total non-perturbative contribution to Ry+A turns out to be 



5NP = (_o.4±1.2) 
5NP = (_o.2±1.2) 



10" 
10" 



(FOPT) , 
(CIPT) . 



(6.7) 



These estimates, despite having errors larger than those quoted previously in the literature, 
must be considered more reliable, as they are the only ones based on an analysis which deals 
explicitly with the theoretical systematic issues discussed above. 

Care must be taken in drawing conclusions from the results of Eq. (6.6). While the 
results do establish that integrated DV contributions to -R^+a are small, it does not follow, 
as repeatedly assumed in the literature, ^^ that DVs can be neglected in the determination 
of as from hadronic r decay data. The reason is that, even if one restricts attention to only 
the quantity Ry+A^ one still needs to determine the D = 6 and D = 8 contributions to S^p. 
This cannot be done in a controlled manner without including values of Sq significantly lower 
than m^ and weights for which integrated DV contributions are certainly not negligible (c/. 
VI and Ref. [29]). Many values of 6p in the literature have been obtained using values 
of 6np taken from analyses with the limitations noted above. In view of the results given 
in Eq. (6.7), the errors on such estimates of 6p are evidently underestimated, often by a 
significant amount. Only improved data will allow these errors to be further reduced. 
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For a recent review, see e.g. Ref. 
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•^M) 



FIG. 5: as(m^) versus x^j covariance matrix reduced by factor 4 (left panel) and factor 9 (right 
panel), V channel, with Smin = 1-5 GeV^. 

VII. FUTURE PERSPECTIVES 

In this section, we speculate on possible improvements relative to the results presented 
in Sec. IV if data with significantly smaller errors were to become available. As mentioned 
before, in principle such data can be extracted from the BaBar and Belle experimental 
results, and it is not unlikely that such an analysis would lead to hadronic spectral functions 
with errors about 2 or 3 times smaller than those of the OPAL data, especially in the upper 
part of the kinematic region, where OPAL statistics are low. 

Therefore, in Figs. 5 and 6, we consider how the results shown in Figs. 1 and 2 would 
change if we used the same central values for the V spectral function as used in Sec. IV, 
but with a covariance matrix scaled by a factor 1/4 or 1/9. We emphasize that this is just 
a speculative exercise. For instance, given that the values of the x^ per degree of freedom 
in Table 1 are of order one, fits with these scaled covariance matrices would give rise to fits 
yielding the same central values, but with poor values of the x^ per degree of freedom. This 
observation refiects just the fact that the data can, of course, not be improved by rescaling 
the covariance matrix, simply because fluctuations in the actual data correspond to the size 
of the actual covariance matrix. 

However, it is of some interest to see what would happen to the properties of the con- 
ditional probability distribution we explored in Sec. IV. The new figures all show that the 
unphysical minimum of Fig. 1 disappears as a function of the rescaling factor, while leaving 
the physical minimum in place. We interpret this as evidence that better data, i.e., data 
with smaller errors, may help resolve the problem that with current data the parameter 5v 
(and therefore, because of the strong correlations, the parameters 71/ and asirnl)) cannot be 
reliably determined without external considerations (c/. Sec. V). Of course, this exercise as- 
sumes that better data would be equally well described by our theoretical parameterization 
of the spectral-function moments we consider in this article. 



VIII. 



CONCLUSION 



In this article, we continued our analysis of hadronic r decays. The main goal is the 
precision determination of the strong coupling at the r mass, asirnl), with good control not 
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n,(m^) 



FIG. 6: Two-dimensional contour plots showing as(m^) versus dy and 7y versus Sy- Top: co- 
variance matrix reduced by factor 4i bottom: covariance matrix reduced by factor 9. Left panel: 
projection onto as{m'^) — 6v plane; right panel: projection onto jy ~ ^V plane. V channel, 
Smin = 1-5 GeV . Blue (darker) areas and green (lighter) areas contain 68%, respectively, 95% of 
the distribution. 



only over statistical, but also over systematic errors. In our previous article, "PI, we presented 
a new framework for such an analysis, in which non-perturbative contributions to the non- 
strange vector and axial hadronic r decays, both from the operator product expansion and 
from violations of quark-hadron duality, can be quantitatively estimated. Since complete 
spectral functions for these channels are available from experiment, the energy dependence 
of these effects can be taken into account. The results of Ref. [27], PI, and the present article 
show unambiguously that it is imperative to take this energy dependence into account in 
order to arrive at a fully consistent understanding of non-perturbative effects. This requires 
the use of a model to parameterize duality violations. We emphasize that a quantitative 
approach to duality violations cannot be avoided, simply because they are clearly present 
in the spectral-function data. The assumption that duality violations are negligible, while 
perhaps reasonable in the past, is no longer acceptable given the current-claimed level of 

• • 22 

precision. 

The specific aim of the present article is two-fold. First, the analysis of VI was based on 



^^ All results reviewed in Ref. [30] claim an error < 0.015 on as{rn^). 
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the 1998 OPAL spectral-function data. The construction of these spectral functions included 
the use of the then available values for the branching fractions for the most important 
exclusive modes. Since these branching fractions are now more precisely known, it is possible 
to update the central values for the spectral functions as well as the corresponding diagonal 
errors. ^^ This update was carried out in Sec. Ill, and the results were subsequently used in 
our fits to the data. 

Our second aim was a more detailed investigation of the quality of the fits that go into 
our analysis. Since these are non-linear, multi-parameter fits, they are of considerable com- 
plexity. The use of a Markov-chain Monte Carlo program made it possible to investigate the 
full posterior probability distribution underlying our most important fits. This allowed us to 
delineate the landscape in parameter space in more detail than through simple minimization. 
This investigation was carried out in Sec. IV. 

Let us summarize what we learned from our new analysis. First, as the reader will note, 
our new results for the value of asi'm'^), contained in Eq. (6.1), are very close to the OPAL 
results of Ref. [2], but with somewhat larger errors. Since the two sets of results correspond 
to sets of data with different normalizations, the near-equality of central values is, in fact, 
purely accidental; as shown in VI, an analysis of the same data as that used by OPAL leads 
instead to significantly smaller values of as(m^). A related observation applies to the errors. 
The errors found in Ref. [2] are smaller simply because systematic effects associated with 
the operator product expansion and duality violations were not considered in that analysis. 
In our opinion, the same observation applies to essentially all determinations of Q;s(m^) from 
hadronic r decays preceding the framework presented in PI. 

Second, while we argued in PI and here that duality violations cannot be reliably left 
out from a quantitative analysis of the spectral functions below the r mass, it turns out that 
the multi-parameter fits thus needed to determine all parameters are at the edge of what 
is possible with currently available data. This is demonstrated in Figs. 1, 2 and 3, which 
show that the probability distributions underlying our fits may have several minima, which 
together span a range of ag values of about 0.27 — 0.34.^"^ Therefore, physical arguments, 
given in Sec. V, are needed in order to narrow down the error on ^^(tti^), and our result (6.1) 
is obtained with the help of these arguments. 

Given this state of affairs, we believe that it would be very interesting to apply our 
analysis to data with much better statistics, which are in principle available from the BaBar 
and Belle experiments. If the non-strange spectral functions that can be extracted from 
these data would be made available, this would allow us to put our analysis framework to a 
much more stringent test. This was demonstrated quantitatively in Sec. VII, where it was 
shown that with much reduced statistical errors one may expect to resolve the ambiguities 
present in the probability distribution constructed from the OPAL data. 

Of course, at present we do not know what the outcome of such an investigation would 
look like. Since fit parameter errors scale as the square root of the scale of the data covariance 
matrix, a factor of 3 improvement in data errors has the potential to produce individual CIPT 
and FOPT fits with errors on ag competitive with those of current lattice determinations. 
Such errors would then be significantly smaller than the difference between current CIPT and 
FOPT results. Theoretical progress on the reliability of various perturbative resummation 
schemes, as embodied in the current discrepancy between CIPT and FOPT, will thus most 



^^ A partial update of the full covariance matrix is also possible, as explained in the Appendix. 
24 This is for FOPT; for CIPT the range is shifted by about 0.02. 
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likely also be necessary. Whether the outcome of such a BaBar- or Belle-based analysis will 
be a more precise determination of the strong coupling near the r mass, or an indication of 
the need to construct more sophisticated representations of non-perturbative effects remains 
to be seen. Either way, we believe that much can be learned from an analysis of the already- 
existing BaBar and Belle data. 
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Appendix A: A partial update of the OPAL spectral functions and covariance ma- 
trices 

OPAL has made publicly available the spectral functions and covariances for the three 
main exclusive modes and inclusive sum over all modes in each of the V and A channels. 
The contributions to the spectral functions corresponding to other exclusive modes (which, 
with the exception of ujtt'tt^, are not measured but constructed using Monte Carlo) are 
not available. The covariances between contributions from different modes are similarly 
unavailable. This limits the extent to which the OPAL inclusive distributions can be updated 
for improvements to the exclusive branching fractions and quantities such as Vud and Be 
which enter the conversion between the inclusive differential decay distributions dBy/A/ds 
and the spectral functions pv/a{,s). 

The procedure for updating Pv/a{.s) was discussed already in the text. The ingredients 
needed for this update are the HFAG branching fractions and the following oo and 77 branching 
fractions, taken from the 2010 PDG compilation: 



5[a;^37r] = 0.892 ± 0.007 , (Al) 

0.0011 

0.0013 ' 



Blu^-K+Ti-] = 0.0153+°-°°ii 



5[?7 ^ vr+TT-vr'^] = 0.2274 ± 0.0028 . 

The latter are needed to convert from the quoted HFAG 37r, 47r and Stt branching fractions 
(corresponding to modes defined such that u and 77 substate contributions are absent) to 
the analogous branching fractions of those exclusive modes tabulated by OPAL (defined 
such that oj and rj substate contributions are included). The corrections to be applied to 
the HFAG branching fractions in order to accomplish this conversion include, in addition to 
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those corresponding to the wrong-current contaminations discussed aheady in the main text, 
those corresponding to the contributions of utt~ to the it' ir^ 7r~ n^ distribution and cjtt^tt^ to 
the 7r~7i~^7c~2n^ distribution produced by the u — )■ 7i'^tt~'it^ decay mode. The remainder of 
the UTT^ contribution represents a mode contribution to be assigned to the V distribution, 
and hkewise, the un^n^ contributions, and the rj-K^-K^ (excluding 7] — )■ tt+tt^tt^) contribution, 
represent mode contributions to be assigned to the A and residual V distributions in the 
OPAL convention, respectively. The remainder of the residual mode contributions consist of 
the wrong-current contamination corrections and (i) for the V channel, the KK, Gtt, KKir 
and KKtttt contributions, and (ii) for the A channel, the 37r~27r"'", 7r~47r°, KKn, KKttti and 
ai (— 7- 7r~7) contributions. We follow OPAL in assuming a fully anti-correlated 50 ± 50% 
breakdown of the KKtt distribution into V and A channel contributions, and employ the 
same assumption for the very small KK2'k contributions not listed by OPAL. The HFAG 
branching fraction for the similarly small ai (— ?■ 7r~7) mode, also not listed by OPAL, has 
also been included in the combined A residual branching fraction sum. 

The inaccessibility of cross-correlations between different exclusive modes limits our abil- 
ity to update the OPAL covariance matrices. We can, however, perform a partial update to 
take into account improvements in the determinations of the constant factors i?e, Sew and 
Vud appearing in the conversion step 

dBv/A{sk)/ds 
Pv/A{Sk) = —^ , (A2) 

where Sk is the midpoint of the A;-th OPAL bin, 

Bk = 12n'^SEw\Vud\^BeWr{yk)/ml = Bwr{yk)/ml , (A3) 

with yk = Sk/ni^ and Wr{y) the (1 + 0) kinematic weight Wr{y) = (1 — i/)^(l + 2|/). From this 
it follows that the relation between the covariances of the spectral function obtained from 
the same dBv/A{s)/ds distribution using new (primed) and old (unprimed) OPAL values for 
the constants Sewi Vud, Be and m^, incorporating also, for completeness, in the updated 
version, the contributions of the uncertainty on rrir neglected by OPAL, is {pi runs over all 
Pv{si) and pA{si)) 



BiBj 

m 




(SpM) + p.p, — - — (A4) 



6m,Y f{-2 + I8yf - 16yf)i-2 + I8y^ - IGyf] 




i^Pi^Pj) 



\mr J V Wr{yi)Wr{yj) ' ^ 

For the current values of the physical quantities appearing in these conversions, we will use 

Sew = 1.0201(3) , (A5) 

iKdl = 0.97425(22) , 

Be = 0.17827(40) , 

TTir = 1.77677(15) GeV , 

from Refs. [35, 36] and [7] for B^ and m,-, respectively. The error on m,- plays no significant 
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role in our analysis. 
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